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ABSTRACT 

We report the detection of a giant planet in a 6.4950 day orbit around the 1.68 M© subgiant 
HD 102956. The planet has a semimajor axis a — 0.081 AU and minimum mass Mp sin i = 0.96 Mj up . 
HD 102956 is the most massive star known to harbor a hot Jupiter, and its planet is only the third 
known to orbit within 0.6 AU of a star more massive than 1.5 M . Based on our sample of 137 
subgiants with M* > 1.45 M© we find that 0.5-2.3% of A-type stars harbor a close-in planet (a < 
0.1 AU) with Mp sin i > 1 Mj up , consistent with hot- Jupiter occurrence for Sun-like stars. Thus, the 
paucity of planets with 0.1 < a < 1.0 AU around intermediate-mass stars may be an exaggerated 
version of the "period valley" that is characteristic of planets around Sun-like stars. 
Subject headings: techniques: radial velocities — stars: individual (HD 102956) — planets and satellites: 
formation 
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1. INTRODUCTION 

The current state of knowledge of planets around 
intermediate- mass (M* > 1.5 M©) stars is reminiscent 
of the general knowledge of exoplanets in 2001. At 
that time there were 32 planets known, mostly orbit- 
ing Sun-like stars. The distributions of semimajor axes 
and masses (Mp sin i ) of these early exoplanet discov- 
eries, drawn from the Exoplanet Orbit DatabasqJ, are 
shown in Figure [T] illustrating the prevalence of giant 
planets in unexpectedly close-in orbitfl As of 2010 May 
there are 31 planets known to orbit intermediate-mass 
stars (M* > 1.5 M©), and the semimajor axes of plan- 
ets in this new class was surprising, but for a different 
reason: there are no planets orbiting c loser than 0.6 AU 
A Johnson et al.ll2007t ISato et alj|2008j ). As iBowler et~aTl 
(2010) showed, the planet populations orbiting host stars 
on either side of 1.5 M© are distinct at the 4-cr level. 

Close-in, Jovian planets are relatively easy to detect 
using radial velocities because of the larger amplitude 
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they induce and the increased number of orbit cycles per 
observing time baseline. Their absence therefore cannot 
be due to an observational bias given the large population 
of planets at longer periods. Instead, it appears that 
stellar mass has a dramatic effect on the semimajor axis 
distribution of planets. However, it is not clear whether 
this effect is a reflection of the process of planet formation 
and migration, or instead related to the effects of the 
evolution of the host stars. 

Stellar evolution may be an important factor because 
Doppler surveys of intermediate-mass stars are largely re- 
stricted to post-main-sequence targets. While massive, 
main-sequence stars are poor Doppler tar gets due to their 
rapid rotation (Vrotsini > 50 km s" 1 : lLagrange et all 
(2009)), their evolved counterparts on the giant and sub- 
giant branches are much slower rotators (U ro tsini < 
5 km s _1 ) and therefore have the narrow absorption 
lines required for precise Doppler measurements. How- 
ever, as stars evolve their atmospheres expand and may 
enc roach upon the o rbits of their planets. Simula tions 
bv INordhaus et al.l (|2010l) , ICarlberg et all (pOOl ) and 
iVillaver &: Livi o (2009) have suggested that the engulf- 
ment planets by the expanding atmospheres of stars can 
account for the lack of clo se-in planets aro und K-giants 
and clump-giants (see also lSato et al J 120081 ) . 

While Doppler surveys have encountered a barren re- 
gion around A stars inward of 0.6 AU, transit surveys 
have discovered two examples of hot Jupiters around 
intermediate-mass stars. OGLE2-TR-L9 and WASP-33 
are 1.5 M© stars orbited by Jovian planets with semi- 
major axes a = 0.041 AU and 0.026 AU resp ectively 
(|Snellen et al.l[200l IColTier Cameron et al.ll2010f ). These 
detections demonstrate that close-in planets exist around 
A-type dwarfs, adding additional concern that the lack 
of planets close to evolved intermediate-mass stars is the 
result of stellar engulf ment. Unfortunately, the com- 
plicated observational and selection biases inherent to 
ground-based, wide-field transit surveys make it difficult 
to measure accurate occurrence rates that can be mean- 
ingful ly compared to th ose measured from Doppler sur- 
veys (jGaudiet al.ll2005l) . 
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Fig. 1. — The semimajor axes and minimum masses (Mp sin i ) of 
the first 32 Doppler-detected planets (gray circles). The majority 
of these planets orbit stars with M* < 1.5 Mq, and the planets 
span a wide range of semimajor axes. Also shown are the first 
31 Doppler-detected planets around massive stars (M* > 1.5 Mq; 
open diamonds), and the sole hot Jupiter, HD 102956 b (filled five- 
point star). Compared to the population of planets around Sun-like 
stars, there is a notable paucity of planets inward of 0.6 AU around 
the intermediate-mass stars. 



Among the various types of evolved stars, subgiants of- 
fer a unique view of the population of hot Jupiters around 
intermediate-mass stars. The radii of subgiants have in- 
flated by only a factor of ~ 2 compare d to their main- 
sequen ce values. The simulations of iVillaver fe Livid 
(2009) show that planets with a > 0.1 AU should be safe 
from the tidal influence of stars near the base of the RGB 
(see their Figure 2). It is only after stars start to ascend 
the RGB and have their radii expand to an appreciable 
fraction of an AU that tidal influences become impor- 
tant. The occurrence rates and semimajor axis distribu- 
tion of close-in planets around subgiants should therefore 
be representative of the properties of planets around A- 
type dwarfs. 

We are conducting a Doppler survey of intermediate- 
mass subgiants at Keck and Lick Observatories to 
study the effects of stellar mass on the physical prop- 
erties and orbital architectures of planetary systems. 
Our survey has resulted in the detection of 14 plan- 
ets around 12 intermediate-mass (M* > 1.5 M Q ) stars, 
and 4 additional planets ar o und less-massive subgiants 
(I Johnson et all I2006L [2007J r200l iBowler et al.l l20Tot 
I Peek et all 1200 9; John son et al.l l2010al Johnson et al. 
2010b, submitted). In this Letter we report the first 
Doppler-detected planet within 0.6 AU of a "retired" 
(former) A star: a hot Jupiter around a 1.68 M Q sub- 
giant. 

2. STELLAR PROPERTIES, RADIAL VELOCITIES 
AND ORBIT 

HD 102956 (=HIP 57820) is listed in the Hipparcos 
Catalog with V = 8.02, B - V = 0.971, a parallax- 
based distance of 126 pc, a nd an absolute magni- 
tude M v = 2.5 (lESAI I1997D . Like most subgiants, 
HD 102956 is chromospherically-quiet with an average 
S = 0.17 ± 0.02 and logfiW = -5.09 on the Mt. 
Wilson scale (|Wright et ai l 12001. We used th e LTE 
spectral synthesis described bv lValenti fc Fischer] (|2005| ) 
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Fig. 2. — Doppler measurements of HD 102956 from Keck Obser- 
vatory phased at the orbital period of the planet. The error bars 
are the quadrature sum of the internal measurement uncertainties 
and 5 m s — 1 of jitter. The dashed line shows the best-fitting orbit 
solution of a single Keplerian orbit. The solution results in residu- 
als with an rms scatter of 6.0 m s _1 and ^/xi = 1-35, indicating 
a good fit to the data. 



and iFischer fc Valentj (|2005fl to estimate the spectro- 
scopic properties. To constrain the low surface gravi- 
ties of the ev o lved s tars we used the iterative scheme of 
iValenti et aLI (J2009D . which ties the SME-derived value 
of log g to the gravity inferred from the Yonsei-Yale 
(Y 2 ) stellar models. Our SME analysis gives T e g — 
5054 ± 44 K, [Fe/H] = +0.19 ± 0.04, log 5 = 3.5 ± 0.06 
and Vrotsini = 0.30 ± 0.5 km s — 1 . We compared the 
star's temperature, luminosity and metallicity to the Y 2 
stellar model grids to estimate a stellar mass M* = 
1.68 ± 0.11 M Q and radius i?* = 4.4 ± 0.07 R Q . All 
of the stellar properties are summarized in Table [TJ 

We began monitoring the radial velocity of HD 102956 
in 2007 April and we have gathered a total of 22 mea- 
surements. After two seasons of observing we noticed 
RV variability with an rms scatter of 41 m s _1 , which is 
much larger than the scatter predicted by the measure- 
ment uncertainties and jitter levels typical of subgiant s 
(|Fischer et al. I [20CJ1 IWrTghtll2005t [Johnson et al.ll2010al) . 
Table [2] lists our RV measurements, ti mes of observa- 
tion and internal errors (without jitter). Uohnson et al.l 
(2010a) estimate a typical jitter of 5 m s _1 based on 
their analysis of 382 RV observations of 72 stable sub- 
giants observed at Keck with HIRES. We add this jitter 
estimate in quadrature to the internal measurement er- 
rors to ensure proper weighting of the data in our orbit 
analysis. 

To search for the best-fitting orbit we used the 
partially-li nearized Keplerian fittin g code RVLIlfl de- 
scribed by iWright fc Howardl ([2009). As an alternative 
to a periodogram analysis, we first stepped through a 
grid of orbital periods sampled from 1 to 100 days in 5000 
equal, logarithmically-spaced intervals. At each step we 
fixed the period, searched for the best-fitting orbit and 
recorded the resulting^ yx^. We found a minimum 

9 http://cxoplancts.org/code/ 

10 We use wxi to indicate the factor by which the observed 
scatter about the best-fitting model differs from our expectation 
based on the measurement errors. Thus, the scatter about our 
model is a factor of 1.35 larger than our average error bar. 



TABLE 1 
Stellar Properties and Orbital Solution for HD 102956 



Parameter 



Value 



V 

B-V 

Distance (pc) 

M v 

[Fe/H] 

T eff (K) 

Vrot sin i (km s 1 ) 

M* (M Q ) 
R* (R©) 
L* (R©) 
logHW 
Age (Gyr) 

P (days) 

K (ms" 1 ) 

e 

Tp (Julian Date) 

u) (degrees) 
Mpsini (M Jup ) 
a (AU) 
R*/a 

Nobs 

rms (ms — 1 ) 



8.02 ± 0.05 

0.971 ±0.01 

126 ± 13 

2.5 ±0.2 

±0.19 ±0.04 

5054 ± 44 

0.30 ±0.5 

3.5 ±0.06 

1.68 ±0.11 

4.4 ±0.1 

11.6 ±0.5 
-5.09 

2.3 ±0.5 

6.4950 ± 0.0004 

73.7 ±1.9 
0.048 ± 0.027 
2455346 ± 0.7 

12 ±40 

0.96 ± 0.05 

0.081 ± 0.002 

0.253 ± 0.008 

22 

6.0 

1.35 



(VxF = 1-15) at periods near 6.5 days. The next lowest 

minimum (yxF — 1-98) occurs at P = 3.47 days. How- 
ever, with a best-fitting eccentricity of 0.95, the orbit 
solution is clearly unphysical. 

We then allowed the period to float in our RVLIN 
analysis with an initial guess of P = 6.5 days, and 
found that a single-planet Keplerian model with a period 
P = 6.4950 ± 0.0004 days, eccentricity e = 0.048 ± 0.027 
and velocity semiamplitude K = 73.7il.9ms~ 1 . The fit 
produces RV residuals with a root-mean-squared (rms) 



scatter of 6.0 m s and reduced 



= 1.35, indicating 



an acceptable fit. W e used the false-alarm analysis of 
Howard et all (12009ft to calculate FAP < 0.001 (see also 
Johnson et al.l I2010b( ) . The resulting minimum planet 



mass is Mp sini = 0.96 Mj up , and the semimajor axis is 
a = 0.081 AU. The best-fitting solution is shown in Fig- 
ure [5J where the plotted error bars are the quadrature 
sum of internal errors and 5 m s _1 of jitter. 

After identifying the best-fitting model, we use a 
Markov-Chain Monte Carlo (MCMC) algori thm to esti- 
mate the parameter uncertainties (See , e.g. iFordl 120051 : 
IWinn et al.l [20071 IBowler et al.l l201(ift . MCMC is a 
Bayesian inference technique that uses the data to ex- 
plore the shape of the likelihood function for each pa- 
rameter of an input model. At each step, one parameter 
is selected at random and altered by drawing a random 
variate from a normal distribution. If the resulting \ 2 
(not reduced by the number of free parameters v) value 
for the new trial orbit is less than the previous \ 2 value, 
then the trial orbital parameters are added to the chain. 
If not, then the probability of adopting the new value is 
set by the difference in \ 2 from the previous and cur- 
rent trial steps. If the current trial is rejected then the 
parameters from the previous step are adopted. We ad- 
justed the width of the normal distributions from which 



the steps are drawn until we achieved a 30-40% accep- 
tance rate in each parameter. The resulting "chains" of 
parameters form the posterior probability distribution, 
from which we select the 15.9 and 84.1 percentile levels 
in the cumulative distributions as the "one-sigma" con- 
fidence limits. In most cases the posterior probability 
distributions were approximately Gaussian. The orbital 
parameters and their uncertainties are listed in Table [1] 

As an additional check on the nature of the RV 
variations, we acquired photometric observations of 
HD 102956 with the T3 0.4 m automatic photometric 
telescope (APT) at Fairborn Observatory. A brief de- 
scription of the T3 data acquisition and reduction pro- 
cedures, as well as the utility of the APT observations 
for eliminating false p ositive detections, can be found in 
I Johnson et"aTT(|2010aft . 

The APT collected two dozen observations in the John- 
son V and B photometric bands between 2010 May 
30 and June 23, near the end of the 2010 observing 
season. The V observations scatter about their mean 
with a standard deviation of 0.0037 mag, consistent with 
T3's measurement precision for a single observation (e.g., 
iHenry et al.l 12000 . Tables 2-3). A least-squares sine fit 
on the 6.4948-day radial velocity period yielded a semi- 
amplitude of 0.0012 ± 0.0010 mag. Identical results were 
obtained for the B observations. This tight upper limit 
to photometric variability provides strong support for the 
planetary interpretation of the radial velocity variations. 

3. THE FRACTION OF INTERMEDIATE-MASS 
STARS WITH HOT JUPITERS 

Our Keck and Lick Doppler surveys of subgiant stars 
have time baselines of 3 and 6 years, respectively, and the 
majority of the target stars have more than 6 observa- 
tions with jitter-limited measurement precision ranging 
from 3-6 m s _1 . Our precision, cadence and time base- 
line provide us with the opportunity to assess the frac- 
tion of intermediate-mass stars with hot Jupiters, which 
we define as planets with a < 0.1 AU. This definition is 
somewhat arbitrary, but it is consistent with definitions 
widely used by other studies of hot Jupiters, which typi- 
cally focus on Solar-mass stars and periods P < 10 days. 
The most comprehens ive study of this kind is that of 
ICumming et al.l (|2008h . who measure an occurrence rate 
of / = 0.004 ± 0.003 for a < 0.1 AU, M P sini > 1 M Jup 
and primarily stars with masses M + < 1.4 M Q . Their 
reported planet occurrence is consistent with / < 0.01 
at 95% confidence. 

Within our sample we restrict our analysis to stars 
with M* > 1.45 M Q (the evolved counterparts of A-type 
stars), iVobs > 3, Vrotsini < 20 km s _1 and no evidence 
of double- lines indicative of an SB2. Our sample contains 
137 stars that meet these criteria. For each star we first 
perform a periodogram analysis and use RVLIN to search 
for orbit solutions near the stro ngest periodicities using 
the same technique described bv lMarcv et al.l (|2005ft . We 
then evaluate the FAP by estimating th e likelihood of 
improving xl over that of a linear fit (see I Howard et al.l 
120091 for further details). For solutions with FAP < 0.01, 
we record the rms of the residuals about the best-fitting 
orbit. For larger FAP values we record the rms about 
the best-fitting linear fit to the RVs. 

Next, we measure the largest velocity semiamplitude 
K up that is consistent with the observed rms scatter for 



TABLE 2 
Radial Velocities for HD 102956 



10.0F 



2.25 



Period [days] 
4.14 6.38 



1.91 



JD 
-2440000 



RV 

(m s- 1 ) 



Uncertainty 
(m s — 1 ) 



14216.805 


11.40 


14429.150 


5.77 


14866.114 


37.10 


14986.851 


-58.52 


15014.772 


59.14 


15016.870 


-49.20 


15284.917 


-79.04 


15285.993 


-14.32 


15313.915 


43.48 


15314.828 


-10.25 


15343.795 


-66.93 


15344.885 


3.98 


15350.786 


-40.57 


15351.880 


34.14 


15372.804 


30.58 


15373.796 


-42.13 


15374.753 


-98.50 


15376.783 


-28.67 


15377.812 


30.78 


15378.749 


50.22 


15379.787 


0.00 


15380.805 


-70.40 



1.18 
1.39 
1.30 
1.27 
1.07 
1.13 
1.34 
1.25 
1.25 
1.18 
1.17 
1.35 
1.27 
1.19 
1.16 
1.22 
1.10 
1.12 
1.10 
1.12 
1.13 
1.13 



simulated plan ets of various pe r iods. Our method is sim- 
ilar to that of iLagrange et al.l (|2009f ) and iBowler et al.l 
(2010). For each star we sample a range of orbital periods 
corresponding to semimajor axes 0.04 < a < 0.10 AUFH 
At each fixed period we generate a sample of 3000 simu- 
lated orbits with random phases and circular orbits sam- 
pled at the actual times of observation. For each sim- 
ulated orbit we also add 5 m s^ 1 of random noise to 
simulate jitter. We then record the distribution of 3000 
simulated velocity rms values and compare the distribu- 
tion to the rms of the measurements. Finally, we adjust 
K until the measured rms is less than that of 99.7% of 
the simulated orbits, and record the semiamplitude as 
K up at that period. 

Repeating this procedure for all 137 stars provides a 
measure the completeness of our survey for planets above 
a given Mp sin i at each semimajor axis sampled in our 
simulations. Figure |3] shows contours of constant com- 
pleteness for our sample, together with the position of 
HD 102956, demonstrating that we are 95% complete for 
M P sini > 1 Mj up and a < 0.1 AU. 

Since each detection or nondetection of a hot Jupiter 
represents a Bernoulli trial, the fraction of stars with 
planets in our sampled range is given by the binomial 
distribution P(f\k,N) ex f k (l ~ f) N - k , where N = 137 
is the number of target stars containing k detections. 
For our sample we measure an occurrence rate of / = 
1.2±o.?, or / < 0.034 at 95% confidence for M P sini > 
0.9 Mj up , which includes a single detection, HD 102956b 
with Mp sini = 0.96 ± 0.05. Restricting our analysis to 
planets more massive than 1 Mj up we find / < 0.025 at 
95% confidence. 

Based on the ir surv ey of main-sequence A-type stars, 
ILagrange et al.1 (|2009t l reported no planets with P < 

11 The lower limit of 0.04 AU corresponds to a/R* = 2 for 
a typical M* = 1.7 Mq subgiant with R* = 4 Rq, which is the 
smallest scaled semimajor axis among the known hot Jupiters listed 
in the Exoplanet Orbit Database. 
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Fig. 3. — Detection limits of the Keck Doppler survey of 
intermediate-mass subgiants as a function of semimajor axis. The 
contours indicate the minimum planet mass detectable with a given 
fractional completeness, which is labeled to the left of each contour. 
The vertical dashed line shows the semimajor axis that is roughly 
equal to 2 stellar radii, assuming R+ = 4 Rq, which is typical of 
stars in the Keck survey. The position of HD 102956 b is marked 
with a solid circle. 

10 days among a sample of 50 startQ, or / < 0.05 
at 95% confidence. However, because their early-type 
stars are such rapid rotators, their achievable velocity 
precision only allowed them to rule out planets with 
Mp sini > 8 Mj up . Thus, our estimate of the fraction of 
massive stars with hot Jupiters represents a significant 
refinement compared to the previous constraints. How- 
ever, our sample size is still too small to provide a mean- 
ingful comparison with the planet fraction measured by 
iCumming et al.l ((2008) for less massive stars. An ex- 
tension of our subgiants survey to the Southern sky is 
therefore warranted, and is currently underway at the 
Anglo Australian Observatory (R. Wittenmyer, private 
communication) . Additional constraints w ill be provided 
by th e Kepler space-based transit survey (jBorucki et al.l 
l200l . 

4. SUMMARY 

We report the discovery of a hot Jupiter (e = 0.048 ± 
0.027, M P sini = 0.96 M Jup , P = 6.4950 days) orbiting 
the subgiant HD 102956. At M* = 1.68 M , HD 102956 
is the most massive star known to harbor a hot Jupiter, 
and this short-period system is the first detected as part 
of a Doppler survey of evolved, intermediate-mass stars. 

The existence of this planet demonstrates that the ob- 
served population of close-in planets around subgiants is 
largely representative of the "primordial" planet popu- 
lation, in that close-in planets have not been adversely 
affected by the relatively mild, post-main-sequence ex- 
pansion of their host stars. This aspect of subgiants, to- 
gether with their jitter-limited RV precision of « 5 m s~ 1 , 
makes them ideally suited for studying the properties 
of short-period planets with a wide range of minimum 



12 Lagrange et al. (2009) do not report stellar masses, but in- 
stead classify stars based on B — V colors. We selected stars from 
their Table 4 with colors 0.0 < B — V < 0.3 as representative of 
main-sequence A-stars in our mass range (Af* > 1.45 Mq). 



masses around intermediate-mass stars. 

Based on our current stellar sample we estimate that 
0.5-2.3% of A-type stars harbor a planet with a < 0.1 AU 
and Mp sini > 1 Mj up , co mpared to the 0-4% occ urrence 
rate around Sun-like stars (jCumming et al.ll2008l ). While 
planets with a < 1 AU are unusually rare around A stars, 
it is possible that there exists a population of hot Jupitcrs 
(a < 0.1 AU) around intermediate-mass stars compara- 
ble to what is found around Sun-like stars. If so, then 
the close-in desert around A stars may simply be an ex- 
aggerated version of the "period valley" observed around 
Sun-like stars, marked by a deficit of planets with periods 
ranging from roughly 10-100 days, a sharp increase in the 
number of detect ed planets beyond 1 AU, and a pile-up 
near P = 3 day dUdrv fc Santos! 120071 iCumming et al.l 
120081 I Wright et al.ll2009f ). Do"p~pler surveys of a larger 
number of massive subgiants, together with careful anal- 
yses of detections from transit surveys, will test this pos- 
sibility. 
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